Polarizabilities of Si 2+ : a benchmark test of theory and experiment 
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We have calculated electric-dipole polarizabilities of the 3s 2 1 5o, 3s3p 3 Po, and 3s3p x Pi states 
of the Si 2+ ion using recently developed configuration interaction + all-order method. Detailed 
evaluation of the uncertainties of the final results is carried out. Our value for the ground state 
electric-dipole polarizability 11.670(13) a.u. is in excellent agreement with the resonant excitation 
Stark ionization spectroscopy value 11.669(9) a.u. [Komara et al, J. Phys. B 38, 87 (2005); Mitroy, 
Phys. Rev. A 78, 052515 (2008)]. This work represents the most precise benchmark test to date of 
theory and experiment in divalent atoms. The near cancellation of the ns 2 1 So ground state and the 
lowest nsnp 
in Si 2+ ion. 



3 Po polarizabilities previously observed in B + , Al + , In + , Tl + , and Pb 2+ is also found 



PACS numbers: 31.15.ac, 31.15.ap, 31.15.am, 06. 30. Ft 



I. INTRODUCTION 

The atomic dipole polarizability describes the first- 
order response of an atom to an applied electric field. 
Atomic polarizabilities have been the subject of consid- 
erable interest and heightened importance in recent years 
due to a number of applications, including development 
of next-generation optical atomic clocks, optical cooling 
and trapping schemes, quantum information with atoms 
and ions, tests of fundamental symmetries, studies of cold 
degenerate gases, thermometry and other macroscopic 
standards, study of long-range interactions, and atomic 
transition rate determinations An imperfect knowl- 
edge of atomic polarizabilities is one of the the largest 
sources of uncertainty in the new generation of optical 
frequency standards 0, [IJ . 

Most of the applications listed above involve monova- 
lent or divalent atoms and ions. There are a number 
of high-precision benchmark tests of experimental and 
theoretical values for the polarizabilities of monovalent 
systems (H. l3T-fl2l] . However, there are few high-precision 
experimental data for the polarizabilities of divalent sys- 
tems which are of particular interest to optical clock de- 
velopment [l3l - fl6| and quantum information [l7j . Most 
recent data for polarizability and Stark shifts of divalent 
systems are compiled in Tables 11, 13, and 14 of Ref. [l|. 

For first row monovalent systems, such as Li and Be + , 
the highest precision determination of the polarizabilities 
by theoretical and experimental methods are found to be 
in good agreement (see recent reviews P, [l|| and refer- 
ences therein). Here, we provide such a comparison for 
a second row divalent species, Si 2+ . We believe that the 
experimental data for this ion provides the most precise 
value of the polarizability of any atomic system with two 
valence electrons. The theoretical calculation presented 
in the present work has the lowest uncertainty reported 
for any divalent polarizability. 



A. Experimental determination of polarizabilities 
from Rydberg spectra 

The polarizability of an ion can be extracted from 
the energies of the non-penetrating Rydberg series of 
the corresponding parent system (see [19j and references 
therein). The polarization interaction between the ionic 
core and the Rydberg electron shifts the energy levels 
away from their hydrogenic values. If the Rydberg elec- 
tron is in a high angular momentum state, it has negli- 
gible overlap with the core. In such cases, the polariza- 
tion interaction provides the dominant contribution to 
the energy shift. This effect is utilized in resonant ex- 
citation Stark ionization spectroscopy (RESIS) (Tfj| - [2c| . 
RESIS experiments have been extremely successful in 
high-precision determination of the ground state polariz 
abilities of H+ and D+ Ne+ fij, Na-like Mg+ " 
Na-like Si 3 + [H, Mg-like Si 2+ fJJ Zn-like Kr 6 + 
Ba+Q, Hgdike Pb 2+ H, Fr " like Th3+ & and Rn " 
like Th 3+ [23|. Quadrupole polarizabilities and transi- 
tion matrix elements have also been determined for some 
of these systems. 

To the best of our knowledge, RESIS experiments pro- 
vide the most precise values known to date of the polar- 
izability of any divalent atomic system. The most pre- 
cise measurement has been carried out for the 3s 2 1 <S'o 
ground state of the Si 2+ ion; a?o = 11.666(4) a.u [20| . 
Later analysis of the RESIS data that included addi- 
tional terms in the polarization expansion yielded cvo = 
11.669(9) a.u [Hj]. Therefore, the Si 2 + RESIS experiment 
presents an excellent opportunity for a high-precision 
benchmark comparison of theory and experiment. 

In this work, we use a recently developed configura- 
tion iteration (CI)+ all-order method |3CH32l| to calcu- 
late properties of Si 2+ . Our value for the ground state 
electric-dipole polarizability of 11.670(13) a.u. is in ex- 
cellent agreement with the RESIS result. Our previous 
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calculation of the Hg-like Pb 2+ ground state polarizabil- 
ity [33| was also in agreement with the RESIS value (ac- 
curate to 0.6%) within our estimated accuracy. 

We note that Mg-like Si 2+ is a particularly interesting 
test system due to its similarity with Mg-like Al + , which 
was used to construct an optical clock with a fractional 
frequency uncertainty of 8.6 x 10 -18 [35j], the smallest 
such uncertainty yet attained. At room temperature, one 
of the largest contributions to the uncertainty budget of 
this clock is the blackbody radiation (BBR) shift. The 
BBR frequency shift of a clock transition is related to 
the difference of the static electric-dipole polarizabilitics 
between the two clock states (36[. We have recently cal- 
culated this effect in Al + using the same Cl+all-order 
approach. Excellent agrement of our present calculation 
with the experiment in the ground state of Si 2+ provides 
an additional test of the approach. At the present time, 
there are no experimental data on the polarizabilities of 
the excited states of Si 2+ to the best of our knowledge. 



II. METHOD 

To evaluate uncertainties of the final results, we carry 
out three calculations in different approximations: CI 
(37|, Cl+many-body perturbation theory (MBPT) [Hj], 
and Cl+all-order [30h32| |. These methods have been de- 
scribed in a number of papers [3(1 H3, HI] and wc 
provide only a brief outline of these approaches and a 
few details relevant to this particular work. 

Our point of departure is a solution of the Dirac-Fock 
(DF) equations 

H ip c = e c ip c , 

where Hq is the relativistic DF Hamiltonian [3l|, HH and 
tp c and e c are single-electron wave functions and ener- 
gies. Self-consistent calculations were performed for the 
[ls 2 2s 2 2p 6 ] closed core, and the 3s, 3p, 3d, 4s, 4p, and 
Ad orbitals were formed in this potential. We constructed 
the B-spline basis set consisting of N = 35 orbitals for 
each of the s, P1/2, P3/2: ■■■ partial waves up to I < 5. 
The basis set is formed in a spherical cavity with radius 
60 a.u. The CI space is effectively complete and includes 
23 orbitals for each partial wave with I = . . . 4. 

The wave functions and the low-lying energy levels are 
determined by solving the multiparticle relativistic equa- 
tion for two valence electrons 371 ]: 



H eS (E n )<f> n = E n $ n . 

The effective Hamiltonian is defined as 

H cS {E) = H FC + £(£), 

where -Hfc is the Hamiltonian in the frozen-core approxi- 
mation. The energy-dependent operator S(-E) takes into 
account virtual core excitations. It is zero in a pure 
CI calculation. The S(-E) part of the effective Hamil- 
tonian is constructed using second-order perturbation 



theory in the CI+MBPT approach [38[ and linearized 
coupled-cluster single-double method in the Cl+all-order 
approach [3l[ . Construction of the effective Hamiltonian 
in the CI+MBPT and Cl+all-order approximations is de- 
scribed in detail in Refs. [3l],[38|]. The dominant part of 
the Breit interaction is included as described in Ref. [3!| . 

The scalar polarizability ao is separated into a valence 
polarizability aJJ, ionic core polarizability a c , and a small 
term a vc that modifies ionic core polarizability due to the 
presence of two valence electrons. The last two terms are 
evaluated in the random-phase approximation (RPA). 
Their uncertainty is determined by comparing the DF 
and RPA values. The small a vc term is calculated by 
adding vc contributions from the individual electrons, i.e. 
a„ c (3s 2 ) = 2a uc (3s), and a vc (3s3p) = a„ c (3s) + a vc (3p). 

The valence part of the polarizability is determined 
by solving the inhomogeneous equation in valence space, 
which is approximated as [4p| ] 

{E v ~ flrfr)|*(v, M')) = D cfi \V (v, J, M)> (1) 

for the state v with total angular momentum J and pro- 
jection M. The wave function ^(v, M') is composed of 
parts that have angular momenta of J' — J, J± 1 that al- 
lows us to determine the scalar and tensor polarizability 
of the state \v, J, M) [4(]]. The effective dipole operator 
D e g includes RPA corrections. 

Unless stated otherwise, we use atomic units (a.u.) for 
all matrix elements and polarizabilities throughout this 
paper: the numerical values of the elementary charge, 
e, the reduced Planck constant, Ti — h/2ir, and the 
electron mass, m e , are set equal to 1. The atomic 
unit for polarizability can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xI0- 8 a (a.u.), where the con- 
version coefficient is Aire^a^/h and the Planck constant 
h is factored out in order to provide direct conversion 
into frequency units; ao is the Bohr radius and eo is the 
electric constant. 



III. RESULTS 

Comparison of the energy levels (in cm -1 ) obtained 
in the CI, CI+MBPT, and Cl+all-order approximations 
with experimental values [HJ is given in Table Q] Corre- 
sponding relative differences of these three calculations 
from experiment are given in the last three columns. 
Two-electron binding energies are given in the first row 
of Table HI energies in other rows are measured from 
the ground state. For a few of the levels, the accu- 
racy of the CI+MBPT calculation is already on the 
order of our expected precision. The accuracy of the 
ground state two-electron binding energy is significantly 
improved in the Cl+all-order calculation in comparison 
with the CI+MBPT one; the CI+MBPT value differs 
from the experiment by —123 cm -1 , while our all-order 
value differs from the experiment by only —7 cm -1 (see 
line one of Table QJ . The inclusion of the all-order core- 
valence correlations significantly improves the differences 
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TABLE I: Comparison of experimental [U and theoretical energy levels in cm~'. Two-electron binding energies are given 
in the first row, energies in other rows are given relative to the ground state. Results of the CI, CI+MBPT, and Cl+all- 
order calculations are given in columns labeled CI, CI+MBPT, and CI+A11. Corresponding relative differences of these three 
calculations with experiment are given in the last three columns in %. 
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between the singlet and triplet states. For example, 
the Cl+all-order value of the 3s3p Pi — 3s3p 3 P\ in- 
terval, 30035 cm" 1 , differs by only 4 cm" 1 from the 
experimental value 30031 cm -1 . The corresponding 
CI+MBPT value, 30120 cm" 1 , differs from the experi- 
ment by 89 cm -1 . As a result, the accuracy of the tran- 
sition energies used in the polarizability calculations im- 
proves in the Cl+all-order approach. 

We separated the effect of the Breit interaction by 
comparing the results of the calculations with and with- 
out the Breit. The Breit contribution to the energies is 
very small, 0.01% or less. However, the inclusion of the 
Breit interaction significantly improves the splittings of 
all triplet states. For example, the 3s3p 3 Pi — 3s3p 3 Po 
and 3s3p 3 p2 — 3s3p 3 Po splittings are 136 cm" 1 and 
413 cm" 1 without Breit, respectively. The values of these 
splittings in our final calculations that include Breit, are 
128 cm" 1 and 389 cm" 1 , in excellent agreement with the 
experimental values, 129 cm" 1 and 390 cm" 1 

We note that the transition energies relevant to the 
calculations of the 3s3p 3 Po polarizabilities are more ac- 
curate than the energies relative to the ground state listed 
in Table |U 

While we do not use the sum-over-state approach in 
the calculation of the polarizabilities, it is useful to estab- 
lish the dominant contributions to the final values. We 
combine our Cl+all-order results for the electric-dipole 
matrix elements and energies according to the sum-over- 



states formula for the valence polarizability [l[ : 



3(2 J + 1) 



7 1. 



E n — E,, 



(2) 



to calculate the contribution of specific transitions. Here, 
J is the total angular momentum of the state v, D is the 
electric-dipole operator, and Ei is the energy of the state 
i. The breakdown of the contributions to the 3s 2 1 S'o, 
3s3p 3 Po, and 3s3p 1 P\ scalar polarizabilities ao of Si 2+ 
in a.u. is given in Table HU Absolute values of the cor- 
responding reduced electric-dipole matrix elements are 
listed in column labeled "D n in a^e. The theoretical 
and experimental [34| transition energies are given in 
columns AP t h and AP cxpt . The remaining valence con- 
tributions are given in rows Other. The contributions 
from the core and vc terms are listed in rows a c and a vc , 
respectively. The dominant contributions to ao listed 
in columns ao(A) and ao(B) are calculated with CI + 
all-order energies and experimental [34[ energies, respec- 
tively. The differences between arj(A) and cvo(B) values 
are small due to excellent agreement of the correspond- 
ing transition energies with experiment. We take ao(B) 
results as final. Our study of the Breit interaction shows 
that it contributes only 0.03-0.07% to the ab initio values 
of polarizabilities. 
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TABLE II: Contributions to the 3s 2 'So, 3s3p 3 Po, and 3s3p 'Pi polarizabilities of Si 2+ in a.u. The dominant contributions 
to the valence polarizabilities are listed separately with the correspon ding absolute values of electric-dipole reduced matrix 
elements given in columns labeled D. The theoretical and experimental [341 ] transition energies are given in columns A_E t h and 
A_B CX pt. The remaining contributions to valence polarizability are given in rows Other. The contributions from the core and 
vc terms are given in rows a c and a vc , respectively. The dominant contributions to a?o listed in columns do (A) and ao(B) are 
calculated with CI + all-order and experimental energies [34|, respectively. 
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-3.794 
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39361 


1.074 
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1.429 




3s3p 'Pi 
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70560 


70680 


1.776 


2.178 


2.181 
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-3s4s 'So 


76185 
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0.996 


0.634 


0.634 
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11.654 
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TABLE III: Summary of the results for the 3s 2 'So, 3s3p 3 P , 
and 3s3p 'Pi polarizabilities of Si 2+ in a.u. and the evalua- 
tion of the uncertainties. First three rows give ab initio results 
for valence polarizabilities calculated in the CI, CI+MBPT, 
and CI+all-order approximations. In the CI+A11 (B) calcu- 
lation, theoretical energies are replaced by the experimental 
values for the dominant contributions. The final results listed 
in row Total Qo are compared with other theory [29(] and ex- 
periment [H [2|. 



Method 
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CI (A) 


11.567 


10.353 


13.040 


CI+MBPT (A) 
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11.519(10) 


10.268(9) 


12.545(15) 
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0.162(9) 


0.162(9) 


0.162(9) 


OLvc 


-0.011(2) 


-0.006(1) 


-0.006(1) 


Total ao 


11.670(13) 


10.425(13) 


12.701(17) 


Theory [29] 


11.688 




12.707 


Theory [41] 


11.75 






Expt. [20] 


11.666(4) 






Expt. [20, 29} a 


11.669(9) 







"This value is a result of revised analysis 12911 of the RESIS exper- 
iment [13]. 



IV. EVALUATION OF THE UNCERTAINTY 
AND CONCLUSION 



There are three contributions to the uncertainties in 
the final polarizability values that arise from the uncer- 
tainties in the valence «q, core cv c , and vc a vc polarizabil- 
ity terms. To evaluate uncertainty in the valence polar- 
izabilities, we compare the results of the CI, CI+MBPT, 
CI+all-order calculations with our final CI+all-order cal- 
culation in which energies in the dominant contributions 
are replaced by their experimental values. The results of 
the last two calculations are given in Table [H] in columns 
ao(A) and ao(B). We summarize the results of all four 
calculations in Table IIIII For consistency, we refer to 
these calculations as CI (A), CI+MBPT (A), CI+A11 (A), 
and CI+A11 (B) since only theoretical energies (in the cor- 
responding approximation) were used in the first three 
calculations. We evaluate the uncertainty of the final re- 
sults in two different ways: (1) as the difference between 
the CI+all-order and CI+MBPT calculations, listed in 
row labeled Diff. (All-MBPT), and (2) as the difference 
between the CI+all-order results with theoretical and ex- 
perimental energies, listed in row labeled Diff. (B)-(A). 
We take the largest of the two uncertainties as the final 
uncertainty in the valence polarizability q;q. The uncer- 
tainty analysis is carried out separately for each state. 
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To evaluate the uncertainty in the a c and a vc con- 
tribution to the polarizability, we calculate these terms 
in both DF and RPA approximations. The DF values 
for the a c and a,; C (3s 2 ) are 0.153 a.u. and -0.0086 a.u., 
respectively. The difference between the RPA and DF re- 
sults is taken to be the uncertainty. Uncertainties of the 
core and valence polarizabilities are added in quadrature 
to obtain uncertainties of the final values. 

The final results listed in row labelled "Total ao" 
are compared with other theory (2i| |4l| and experi- 
ment [lo, H?|. Our value for the ground state polariz- 
ability is in excellent agreement with both original RE- 
SIS value [2(| value and the revised RESIS analysis [2!| . 
Our values for the ground and 3s3p 1 P\ state polarizabil- 
ities are in excellent agreement with theoretical values 
obtained with large-scale CI calculation with semiempir- 
ical inclusion of the core polarization 29] . The CI result 
of 41] is consistent with other values; the small difference 
is probably due to omission of the highly-excited states 
in the valence CI and restricted treatment of the core 
excitations in |4l| . 

We note that the values of the 1 5'o and 3 Pq polarizabil- 
ities given in Table IIIII are very similar, their difference 
is only 10% of the ground state polarizability. 



In summary, we have carried out a benchmark test 
of the theoretical and experimental determination of the 
ground state polarizability of the Si 2+ ion. Our final re- 
sult is in excellent agreement with the RESIS experimen- 
tal value [20I . [29( 1 . High-precision recommended values 
are provided for the excited state 3s3p 3 Pq and 3s3p 1 P\ 
polarizabilities. The near cancellation of the ns 2 Sq 
ground state and the lowest nsnp 3 Pq polarizabilities re- 
ported for B + , Al + , In + , Tl + , and Pb 2+ is also observed 
for the Si 2+ ion. 
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